Submitted for publication in The Astrophysical Journal 

Preprint typeset using FiTpX style emulateapj v. 5/2/1 1 



WHAT CAN WE LEARN FROM THE RISING LIGHTCURVES OF RADIO ACTIVELY-POWERED SUPERNOVAE? 

Anthony L. Piro 1 and Ehud Nakar 2 

Submitted for publication in The Astrophysical Journal 

ABSTRACT 

The lightcurve of the explosion of a star with a radius < 10— lOO/?0 is powered mostly by radioactive decay. 
Observationally such events are dominated by hydrogen deficient progenitors and classified as Type I super- 
novae (SNe I), i.e., white dwarf thermonuclear explosions (Type la), and core collapses of hydrogen-stripped 
massive stars (Type Ib/c). Current transient surveys are finding SNe I in increasing numbers and at earlier 
times, allowing their early emission to be studied in unprecedented detail. Motivated by these developments, 
we summarize the physics that produces their rising lightcurves and discuss ways in which observations can 
be utilized to study these exploding stars. The early radioactive-powered lightcurves probe the shallowest de- 
posits of 56 Ni. If the amount of 56 Ni mixing in the outermost layers of the star can be deduced, then it places 
important constraints on the progenitor and properties of the explosive burning. In practice, we find that it 
is difficult to determine the level of mixing because it is hard to disentangle whether the explosion occurred 
recently and one is seeing radioactive heating near the surface or whether the explosion began in the past and 
the radioactive heating is deeper in the ejecta. In the latter case there is a "dark phase" between the moment of 
explosion and the first observed light emitted once the shallowest layers of 56 Ni are exposed. Because of this, 
simply extrapolating a lightcurve from radioactive heating back in time is not a reliable method for estimating 
the explosion time. The best solution is to directly identify the moment of explosion, either through observing 
shock breakout (in X-ray/UV) or the cooling of the shock-heated surface (in UV/optical), so that the depth 
being probed by the rising lightcurve is known. However, since this is typically not available, we identify 
and discuss a number of other diagnostics that are helpful for deciphering how recently an explosion occurred. 
As an example, we apply these arguments to the recent SN Ic PTF lOvgv. We demonstrate that just a single 
measurement of the photospheric velocity and temperature during the rise places interesting constraints on its 
explosion time, radius, and level of 56 Ni mixing. 

Subject headings: hydrodynamics — shock waves — supernovae: general 



1. INTRODUCTION 

A typical supernova (SN) lightcurve is powered by a com- 
bination of two sources: (1) the energy deposited by the SN 
shock and (2) the radioactive decay of 56 Ni that was synthe- 
sized during the explosioifl The first light from shock break- 
out and the following early emission (for the first minutes 
to days) is dominated by the shock-deposited energy. The 
very late emission, after the radiation can efficiently diffuse 
through the entire ejecta (> 100 d), is dominated by radioac- 
tive decay. The relative influence of these power sources dur- 
ing the time in between depends on two main factors: (1) the 
progenitor radius and (2) the amount of 56 Ni. The reason 
for this is that the bolometric luminosity from shock heating 
increases linearly with the progenitor radius (as we discuss in 
82.21 ). while the peak of the radioactively-po wered lumino sity 
is roughly linear with the total mass of 56 Ni (!Arnettlll979l) . 

Core-collapse SNe typically synthesize ~ 0.01 -0.1M Q of 
56 Ni, and the progenitor radii range between ~ 1 - lOOO/?©. 
Therefore, in more extended progenitors the main SN event 
is dominated by shock heating, while in more compact ones 
it is dominated by radioactive power. The most common SNe 
II-P are explosions of red supergiants with ~ 500/?©. They 
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exhibit an extended plateau phase for about a hundred days, 
which is powered by the cooling of shock-heated material. 
The progenitors of the ra re 1987A-like SNe II are blue super - 
giants with R* ~ 50R Q (IWoosleylll98l iKleiser et al.ll201 lh . 
and their lightcurves are dominated by radioactive decay start- 
ing a week after the explosion. Similarly, SNe l ib with com- 
pact p rogenitors (R* < 10 R Q ) like SN 201 ldh dArcavi et alj 
201 1) show evidence of both shock-heated material and a sep- 
arate radioactive peak. The progenitors of SNe Ib/c are mas- 
sive stars that were stripped of their hydrogen envelope and 
are also fairly compact with R* < IQRq. Indeed SNe Ib/c are 
dominated by radioactive power starting a few days (or ear- 
lier) after the explosion. Finally, for the compact white dwarf 
progenitors of SNe la with R* < O.OIRq, the shock-heating 
lightcurves are so dim that they have never been observed, 
and our knowledge of these events is only possible due to the 
synthesis of ~ 0.5M© of 56 Ni. 

Since the physics that governs the lightcurves of most SNe 
with < 5Orv0 has many similarities, a single theoretical 
framework should roughly describe their main qualitative fea- 
tures. Their emission during the first few days is especially 
exciting because it probes the shallowest layers of the progen- 
itor. This can teach us about the exploding star's radius, and 
constrain the surface composition and velocity/density gradi- 
ents that reflect details of the explosive burning. The shock- 
heating contribution in SNe I has been w ell-studied in the lit- 
erature, both with semi-an alytic models dNakar & Sarill20 10b 
iRabinak & Waxmanl201 lh and w ith detailed radiative transfer 
simulations (iDessart et al.11201 ll) . The rising lightcurve from 
radioactive heating has been explored for the cases of SNe la 
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dPiro 12012b and to stud y the impact of 56 Ni mixing in SNe 
Ib/c JPessart et al.ll2012l) . 

Concurrent with these theoretical studies, transient sur- 
veys li ke with the Katzman A utomatic Imaging Telescope 
(KAIT; lFilippenkoeTal]|2001l) and by the Palomar Tran- 
sient Factory (PTF; iRau et al] I2009t lLaw et al.1 12009t) and 
the Panoramic Surv ey Telescope and R apid Response Sys- 
tem (Pan-STARRS: iKaiser etall 120021) are finding increas- 
ing numbers of these events, especially at early times. Best 
known amon g these is SN 20 lif e, the closest SNe la in the 
last 25 years iNugent et alJlMTlt iBloom et al.ll2012l) . Other 
SNe la repo rted within the la st year with e arly data include 
SN 2009ig ( lFolevetal.l l 2"0T2h. SN 2010jn (iHachinger et al.l 
I20T1) . and SN 2012cg Isiiverman et alJ 120121) . SNe Ib/c 
have also been increasing ly well-studi e d in th e optical at early 
times, as summarized by iDrout et aD (1201 lh. Other particu- 
lar recent events in clude SN 2008D (ISoderberg et al.l l2008h 
iModjazetal] 120091) and PTF 12gzk (Ben- Ami et all 120121) . 
Motivated by these exciting developments, we ask: what can 
and what cannot be learned from the early optical lightcurve 
of radioactively-powered SNe ? Being more abundant, we fo- 
cus on parameters that are typical to SNe I, although most of 
our conclusions can also be applied to subclasses of SNe II 
that are dominated by radioactive heating. 

In our study we first consider instances where only a pho- 
tometric lightcurve is available, as is often the case at early 
times. One of the main issues is how much we can infer about 
56 Ni mixing in the outermost layers of the star with this lim- 
ited information. Our main conclusion is that it is hard to de- 
termine the mixing of 56 Ni based on the lightcurve alone with- 
out additional information. The reason is that the lightcurve 
can provide an estimate for the total mass of 56 Ni that is ex- 
posed to the observer at any given time, but it cannot provide 
a good handle on the fractional mass of 56 Ni with respect to 
the total exposed mass (which can be determined if the time of 
explosion is well constrained, e.g., by detection of the shock 
breakout or the cooling envelope phase). Putting it differently, 
there is a degeneracy in the lightcurve between a recent explo- 
sion with a high fraction of 56 Ni in the outermost layers and 
an older explosion where 56 Ni resides only in deep material. 
In the latter case the SN has a "dark phase" that can persist for 
up to a few days after the explosion. During this time no 56 Ni 
is exposed and its shock-heated cooling emission is often too 
faint for detection even when deep observations exist. For this 
reason, a simple extrapolation of the lightcurve to early times 
cannot reliably constrain the time of explosion. 

We then consider what additional information can be ex- 
tracted by using limited spectral data that provide the pho- 
tosphere temperature and velocity. We show that if detailed 
color evolution, or if one or even better two spectra (pos- 
sibly of low signal to noise) are available during the rise, 
then the degeneracy between explosion time and 56 Ni depth 
can be at least partially alleviated. Key among our results 
is equation ( fT9b . which allows one to estimate a lower limit 
on the time of explosion using merely a single simultaneous 
measurement of the bolometric luminosity, temperature, and 
photospheric velocity. The methods we discuss are comple- 
mentary to detailed spectral studies that probe element mixing 
based on their absorption and emission features in the rare 
cases where early high signal to noise sp ectra are available 
(e.g. JSauer et alJl2006t|Parrent et alj|2012l) . 

In the following we begin in $2]by discussing each of the in- 
gredients that shape the electromagnetic emission starting less 
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FIG. 1. — Schematic diagram showing the early lightcurves for SNe I and 
how they relate to the photospheric radius and velocity. The top and middle 
panels demonstrate how the relative positions of the shock-heating lightcurve 
(blue curves), diffusive tail lightcurve (purple curves), and the 56 Ni lightcurve 
(red curves) can differ depending on the depth of the 56 Ni. The total observed 
lightcurve is the sum of these three components. Note that when the 56 Ni 
is deposited deeply (in the middle panel) and the shock-heating lightcurve 
(blue curve) is below the detection limits, there can be a significant dark 
phase between the time of explosion and the moment of first detection. In 
the bottom panel we show how the photospheric radius (orange curve) and 
velocity (green curve) evolve with time. Depending on the position of the 
56 Ni lightcurve, different photospheric radii, velocities, and velocity gradi- 
ents will be present when the rising lightcurve is observed. These provide 
clues about the depth of the 56 Ni and the time of explosion, as summarized 
in the middle panel and discussed in the text. 



than an hour and up to weeks after explosion. This discussion 
provides a useful general guide for interpreting early observa- 
tions of radioactively-dominated SNe. In $3]we investigate a 
specific SN Ic in some detail (PTF lOvgv) as a test case for 
applying these arguments and techniques. We conclude in §|4] 
with a discussion of our results and a summary of important 
conclusions that should help facilitate better constraints from 
future SNe I observations. 

2. EARLY EMISSION FROM TYPE I SUPERNOVAE 

When a hydrogen-poor star explodes as a SN, a few impor- 
tant events occur in the moments before and after the first op- 
tical emission is seen. In this section we summarize the main 
properties of each of these events, their observational conse- 
quences, and how detections of some or all of these events can 
be used to put constraints on the properties of the exploding 
star. To guide the discussion we will be referring to the dia- 
gram in FigureQ] which shows the time-dependent luminosity 
components, photospheric radius, and velocity. 
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2.1 . Shock Breakout 

Just prior to emission, a shock is traveling through the en- 
velope. This heats and accelerates the material, unbinding 
it from the star. A radiation-dominated shock accelerates in 
the decreasing density of the outer edge of the star with a 
shock vel ocity that scale s with the density as v s oc p - ' 3 where 
(3 ps 0.19 (ISakurailll960l) . The shock continues to shallower 
regions until the optical depth falls Iotkc/ v s , where c is the 
light speed. At this point the photons are no longer trapped; 
they stream away and the shock dies. This is what is typically 
referred to as the "shock breakout" UV/X-ray flash, and it has 
been frequently studied because of its strong dependence on 
the radius (which determines both the energy budget of the 
shock and the timescale of the emission), all owing the pro- 
genitor star to be studied from its detection dColgatelll974t 
lFi5H!97aiKlein & Chevalieri[l978l) . 

Recently it was realized that the radiation behind the shock 
falls out of thermal equilibrium for small radii hydrogen- 
stripped massive stars with h igh shock veloc i ties, and that 
the br eakout will be in X-rays dKatz et al.ll2010tlNakar & Saril 
120101) . In SNe la, the shock achieves relativ istic velocity, 
and it s breakout emission is in the MeV range dNakar & Saril 
l2012h . In any case, although shock breakout is the first indi- 
cation of the explosion, its impact for SNe I in the optical/UV 
bands is negligible. For this reason it is not discussed further 
in this paper and not plotted in Figure Q] 

2.2. Shock-Heated Cooling Lightcurve 

Immediately after shock breakout, the observed radia- 
tion is out of thermal equilibrium until rou ghly the time 
when the gas doubles its radius at t ps R*/vf dNakar & Saril 
120101) , where Vf ps 2v t is the final velocity of the ma- 
terial (iMatzner & McKeel 119991) . which is within minutes 
or less. The observed photons then gain thermal equi- 
librium _and_the expmsioji_e^ 

phase dChevalierlll992t iPiro et al.l l20lot iNakar & Saril 120101: 
iRabinak & Waxmanl 1201 ll) . As the material that has been 
heated by the shock expands, a thermal diffusion wave begins 
backing its way through the ejecta. Above the depth of the dif- 
fusion wave, material cools via photon diffusion. Below this 
depth, material evolves adiabatically. For each fluid element, 
there is then a competition between adiabatic cooling and dif- 
fusive cooling that controls the energy density of that material 
at the moment its photons begin streaming out of the star. This 
determines the observational signature of shock-heated coo l- 
ing, leading to a bolometric luminosity (Nakar & Sari 2010IP 1 
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(1) 



(2) 



where E = £51 10 51 erg is the explosion energy, M = M\Mq 
is the ejecta mass, R t = R*/Rq, and k is the opacity with 
kq.2 = k/0.2 cm 2 g" 1 being the canonical value during the 
cooling phase for fully ionized hydrogen-free gas. For all the 
scalings in this paper we assume a polytropic index of « = 3, 
as is relevant for compact, radiative stars forming SNe Ib/c, 



4 The luminosity pie-factor of Nakar & Sari 12010) is divided here by a 
factor of 2.5 following the numerical result of Katz et al. (2012) 



or relativistic, degenerate WDs exploding as SNe la. The key 
result is that the luminosity during the shock-heated cooling 
is directly proportional to the progenitor radius. 

As long as the opacity is dominated by scattering the ob- 
served temperature is determined at the thermalization depth, 
below the photosphere, where the optical depth is r > 1. 
The thermalization optical depth, r c , and the absorption op- 
tical depth, T a b s , satisfy at the thermalization depth T e r a b s ~ 1. 
Equation © assumes that r a b s is dominated by free-free and 
minor correction are expected when b ound-free absorptio n 
dominates (for a detailed discussion, see lNakar & Sarill2012h . 

In the optical/UV, the luminosity is rising as long it is in the 
Rayleigh-Jeans tail, resulting in (see Appendix lAl 



^opt/UV r 



1.5 
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When optical/UV luminosi ty is observed to rise more steeply 
than this (as in PTF 12gzk: lBen-Ami et al.ll2012l) . it indicates 
that the rise from shock-heated cooling is not being seen, and 
the time of explosion is not confidently constrained. This 
phase continues for several hours in a core-collapse SN un- 
til T c crosses the UV (at which point the UV flux starts to 
drop) and reaches the optical band. This is also the point that 
recombination starts in cases that radioactive decay does not 
play an important role. In SNe la this phase is terminated af- 
ter about an hour when the diffusion wave reaches material 
where the shock was not radiation do minated and the cooli ng 
envelope emission drops significantly (IRabinak et al.ll2.0fTh . 

Once T c drops sufficiently in a core-collapse SN, a re- 
combination wave begins backing its way in from the sur- 
face. The temperature drop becomes more gradual, settling at 
ps 5000-8000 K an d the luminosity drop stops or even gently 
rises (as discussed in lDessart et alj[20TTt also see Goldfriend, 
Sari & Nakar, in preparation). As far as optical photome- 
try is considered, both the temperature and the luminosity are 
roughly constant and the SN enters a "plateau phase." This 
plateau is similar to the o ne observed in SN II-P dPopovl 19931 : 
iKasen & Woosle v 2009), except that for SNe I it is dimmer 
and that it starts earlier, both due to the smaller progenitor 
and enhanced adiabatic losses. Setting T r ps 0.6 eV as foun d 
for the plateau in the SN Ib/c models of lDessart et all d2011h . 
equation (f2]) implies that the optical plateau starts at 
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and together with equation (Q]i we find that 



pO.85pO.78 
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ergs 



(4) 



(5) 



is roughly the plateau luminosity 

To conclude, the detection of the shock-heated cooling 
phase in core-collapse SNe I is very challenging given its 
low luminosity and short duration, but the rewards are also 
high. Photometry alone of the rising phase provides tight con- 
straints on the explosion time, and both the optical rising and 
subsequent plateau phase provide constraints on the progeni- 
tor radius (if not overcast by radioactive heating as we discuss 
next). A convenient method to identify that a rising optical 
emission is due to shock-heated cooling emission (and not 

5 Applying the same line of arguments to a red supergiant with R* = 
5OOi?0, M = 15-Mq, and k = 0.34 cm 2 g -1 , results in t p ~ 10 days and 
L p ss 10 42 ergs -1 , compatible with observations of SNe II-P. 
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radioactively powered emission), is that the optical emission 
rises (as predicted by eq. 13), while the bolometric luminos- 
ity drops. This can be seen, for example, by having simultane- 
ous UV and optical coverage in the time that T c sits between 
these two frequency windows. 

2.3. 56 M Shallower than the Diffusion Depth 

If there was no radioactive energy input, this would be the 
end of the story for the electromagnetic signal from SNe, and 
most SNe I would be too dim to have ever been detected. But 
eventually the radioactive decay of 56 Ni starts heating the ex- 
panding ejecta, and this is where their lightcurves begin rising 
in earnest. 

When the thermal diffusion wave reaches the shallowest de- 
posits of 56 Ni, the energy generation from 56 Ni roughly goes 
directly into the observed bolometric luminosity (as shown 
with the red curves in Figure[T|i, so that 



L 56 ps M 56 e, 



(6) 



where is the mass of 56 Ni that is exposed by the diffusion 
wave, and the specific heating rate from 56 Ni is 



e(f) = e Nie -'/' N ' + eco(e 



-'/'Co . 



"f/'Ni 



(7) 



where e Ni = 3.9 x 10 10 erg g" 1 s" 1 , f Ni = 8.76 days, e Co = 
7.0 x 10 9 erg g" 1 s" 1 , and fco = 111.5 days. It should be noted 
that in our more detailed analysis in Appendix [B] we show 
that if L56 is the observed luminosity, then equation (O is 
only accurate up to a factor of ~ 2. Even with this small un- 
certainty, the rough level of 56 Ni mixing that can be deduced 
from the observations is still robust as we will discuss. 

During this phase the depth of the diffusion wave is impor- 
tant because it tells us which part of the exploding star is being 
probed by the observations. The depth in mass is related to the 
time after explosion by (see AppendixO 
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c 51 
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0.32 



1 day 
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where ko.i = k/0.1 cm 2 g -1 is the canonical value we use for 
the radioactively powered phase, since the gas is partially ion- 
ized. For the estimates in this paper, we use quantities that 
roughly correspond to the SNe Ib/c, but the same general ar- 
guments apply to SNe la (as well as the rare SNe II that are 
radioactively powered) with different prefactors. To see how 
much these can vary for other progenitors see Appendixlcl 
Also note that equation (0 becomes less accurate closer to 
the peak of the lightcurve, when AMdiff grows more slowly 
with time. 

When powered by radioactive heating, the observed lumi- 
nosity is no longer sensitive to the progenitor radius as can be 
seen by the lack of an explicit dependence on R* in equation 
©. Instead, it provides a direct measurement of the amount 
of 56 Ni (via eq. J6)) at the location of the diffusion wave. 
If the explosion time is known, equation (J8j provides AMdiff 
into which this 56 Ni is mixed. From this the mass fraction 
X56 ps M56 / AMdiff can be inferred as a function of the depth 
AMdiff. Thus, in principle, detection of the rise of the 56 Ni 
lightcurve should provide an estimate of the distribution of 
56 Ni (as attempted for the SN la 201 lfe bv lPiroll2012l) . which 
is helpful for understanding the nature of the explosive burn- 
ing and the outer structure of the progenitor. 



2.4. 56 M Deeper than the Diffusion Depth: "Diffusive Tail" 

In practice, the exercise described above is not as straight 
forward as one might think. The reason is that if there is a 
steep increase in the 56 Ni abundance, then the assumption that 
the observed emission is generated only by the composition at 
the location of the diffusion wave may not be valid. Instead 
a "diffusive tail" of the energy released in 56 Ni-rich layers 
deeper than the diffusion depth (shown as purple curves in 
Figure [U may actually dominate over the energy released in 
56 Ni-poor layers shallower than the diffusion depth (shown as 
red curves in Figure[T|i. Next we calculate the contribution of 
this diffusive tail. 

Consider a deposit of 56 Ni at some depth d in the star. The 
56 Ni decays to produce gamma-rays, which are absorbed and 
create thermal photons. The photons then spread due to diffu- 
sion, creating roughly a Gaussian distribution around the 56 Ni. 
At each time f , the Gaussian has a width \f~Kt , where K is the 
diffusion coefficient, which is related to the diffusion time by 
fdiff = d 2 /K. For such a distribution, there is a small, but non- 
zero, fraction of photons that escape from the star because 
they have diffused by a distance > d from the 56 Ni depth. The 
fraction of photons that reach the surface at time t is 



Escaping fraction « erfcy7diff72f, (9) 

where erfc is the complementary error function. Therefore, if 
at some time t ' the diffusion wave reaches a layer of 56 Ni that 
would be producing a luminosity L^it 1 ) according to equation 
(O, this implies for previous times t < t ' that the diffusive tail 
from this layer also produces a luminosity. Since ?diff oc pr 2 oc 
1 jt, this diffusive tail luminosity is 



Lm(t<t')=Ls 6 (t') 



e(t) erfc(*7V2*) 
W) erfc(l/V2) 



(10) 



(11) 



A simpler approximation to this scaling is 

, e{t) t e~'' 2/2 ' 2 
Ltoa(t <t )psL 56 (j )-7-pr—. — 1775-, 
e(t') V e y l- 

which is accurate within ps 20% for t'/2 <t <t' even though 
Ltaii changes by a factor of ps 9 over this range of times. 

From this discussion one can see what we meant earlier in 
this section when we mentioned "a steep increase in the abun- 
dance of 56 Ni." Consider the two 56 Ni distributions shown in 
Figure [2] In the top panel the 56 Ni has a rather shallow dis- 
tribution which produces the L^{t ) (shown by the red curve). 
Therefore if we consider the luminosity L56 at some time t', 
and then trace back the diffusive tail implied from that depth 
using equation dTOb (shown as a purple curve), the diffusive 
tail always falls below the lightcurve. This is a case where 
the shallow 56 Ni prevents the diffusive tail from having a no- 
ticeable impact. In the bottom panel the 56 Ni has a steeper 
distribution. Now, when the diffusive tail is drawn back from 
a point at time ?', it exceeds the lightcurve. In this latter 
case the diffusive tail will dominate the observed rise. 

Therefore to evaluate whether the diffusive tail is important 
at a given depth the key quantity to consider is the slope of 
Ltaj]. Using equation ( fill ), the logarithmic derivative is 



d\nL 



tail 



dint 



1+, 7 



(12) 



For t <C t', this shows that L t aii(0 can become rather steep, 
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FIG. 2. — Schematic diagram showing how the slope of the 56 Ni distribution 
determines the relative importance of the diffusive tail. In both the top and 
bottom panels, a luminosity L^{t ') would be produced at time t' from the 
diffusion wave probing the 56 Ni distribution. In the top panel, the slope of 
the 56 Ni is shallow, which produces a shallow slope for L56 (shown by the red 
curve). If a diffusive tail luminosity is drawn back from this point (shown by 
the purple curve), it always falls below the 56 Ni lightcurve. Thus the diffusive 
tail from this depth is not important. In the bottom panel, the slope of the 
56 Ni distribution is steeper. Therefore, the diffusive tail from larger depths 
overpowers the heating from 56 Ni at shallower depths, and the diffusive tail 
from this depth impacts the lightcurve. The logarithmic slope of the diffusive 
tail near the location t ~ t' is roughly 1 + (r'/0 2 ~ 2. 

as should be expected because of the exponential dependence 
it has on time. For t w ?', the slope of L ta ii(f) is at its shal- 
lowest with roughly L ta ii oc t 2 . From this we conclude that a 
rising lightcurve that is shallower than ~ t 2 indicates that di- 
rect heating from 56 Ni is dominating, as shown in the upper 
panel of Figure |2] In such a case, M56 can be approximated 
from the observations and X56 ~ M56 / AMjiff can be inferred 
as a function of time if the explosion time is well constrained. 
Otherwise, a rising lightcurve steeper than ~ t 2 indicates that 
the diffusive tail is having an impact as shown in the bottom 
panel of Figure|2] and this diffusive tail must be accounted for 
in order to derive Xs^. 

This conclusion is only approximate since we are focusing 
on the the diffusive tail from a single depth. This is likely not 
too bad of a simplification, since the diffusive tails from 56 Ni 
at even larger depths are exponentially suppressed. In detail 
though, the luminosity of the diffusive tail should depend on 
the sum of contributions from all depths. We calculate and 
discuss the impact of this in Appendix |B] The conclusion is 
that even in instances when local radioactive heating domi- 
nates, the expected luminosity given by equation © is never 
more accurate than a factor of ~ 2. 

2.5. The Importance of the Time of Explosion 

Our discussion thus far makes it clear that it is very impor- 
tant to know the time of explosion. It sets the time at which 



the thermal diffusion wave begins backing its way into the 
expanding ejecta, and from this at all later times it is roughly 
known what depths of the exploding star are being probed by 
the observations via equation dH). 

When the time of explosion is not known from direct de- 
tection of shock breakout or shock heating of the surface lay- 
ers, our discussion of the early lightcurve should also pro- 
vide some reason for caution. To illustrate the problem, in the 
top and middle panels of FigureQ]we compare the lightcurves 
for different 56 Ni depositions. In the top panel, the 56 Ni is 
deposited into rather shallow layers, therefore the timescale 
between the beginning of the explosion and the rise of the 
lightcurve is fairly short. In this case, the time of explosion 
could be reasonably well-approximated by extrapolating the 
56 Ni lightcurve back in time. In the middle panel, the 56 Ni 
is deposited in deeper layers, and correspondingly the delay 
between the shock heating and rising lightcurve is longer. In 
this case, extrapolating the 56 Ni lightcurve back in time would 
provide a poor estimate for the time of explosion. We are 
lead to the following important conclusion: one cannot sim- 
ply estimate the time of explosion by extrapolating the rising 
lightcurve back in time because its position relative to the mo- 
ment of explosion depends on the depth of radioactive heat- 
ing. This means that for events like SN 201 lfe, t he constraints 
on cannot be as tight as previously reported dBloom et al.l 
2012) when the explosion time is not known. 

The earliest detection of the rising lightcurve from 56 Ni 
does not probe the shallowest layers of the star, but merely 
the shallowest deposits of radioactive heating. For example 
if the diffusive tail only reaches 56 Ni after traveling ~ 0.1M Q 
below the exploding star's surface, then it will take ~ 2 days 
to detect this depth (using eq. JH]). Depending on the specific 
parameters for a given SN, this timescale can even be a few 
times larger. Similar delays between the moment of explo- 
si on and the rising ligh tcurve are seen in the numerical work 
of lDessartetaflfeOllI) . 

2.6. Clues about the Depth of 56 Ni 

When the shock breakout and shock heating are not de- 
tected, additional information from color, or even better spec- 
troscopic, observations can be used to break the degeneracy 
between the depth of 56 Ni and the time of explosion. We sum- 
marize some of the properties of SNe that are most useful for 
doing this in this section and the next. 

In the bottom panel of FigureQ]we schematically show the 
time-dependent radius of the photosphere (orange curve) dur- 
ing a SN . For a polytropic index of n = 3, the photospheric 
radius is (Rab inak & Waxrrianll201 ll) . 
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where we have suppressed the dependence on the density 
structure factor f p and set it to 0.01 (see Appendix ID1. Even 
though r p h is moving into the star as the ejecta expands, it is 
always moving out in an Eulerian frame. The observed color 
temperature is 



T c *( Lt < 
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where r c is the color radius and r c is the optical depth at the 
color (thermalization) depth. The inequality comes from the 
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fact that t c > 1 and r c < r p h. For typical parameters r c is not 
much larger than unity (see Appendix[E]i and so is r p h/r c , so 
this inequality is also a fair approximation of T c . Although 
r p h roughly always scales as a power-law with time (eq. | [T3l ). 
L{t) can vary greatly depending on the 56 Ni deposition. Thus, 
the color evolution during the rising phase provides the first 
clue on the 56 Ni depth. For deep 56 Ni, the initial rising phase 
of L is exponential, so from equation dT4b we see that T c 
should be rising as well. On the other hand, if 56 Ni is shal- 
low then L rises more gradually and T c is roughly constant or 
even slowly decreasing during the lightcurve rise. 

In the bottom panel of Figure [T] we also show the photo- 
spheric velocity (green curve). Taking v p h = r^/t, we find 
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M 1 



0.28 



1 day 



-0.22 



km s 



(15) 



Unlike the photospheric radius, the photospheric velocity is 
decreasing with time. This is because the photosphere is mov- 
ing into the star in a Lagrangian sense, and reaching ma- 
terial that is moving at successively lower velocities. Fur- 
thermore, the time derivative of the photospheric velocity is 
dvph/dt oc f~ L22 . Therefore we expect the change in the pho- 
tospheric velocity to become more gradual with time. 

Using the above arguments, we conclude that there are four 
main ways in which the depth of 56 Ni heating may be qual- 
itatively inferred. All other things being equal, a larger 56 Ni 
depth implies the following. 

1 . The photospheric radius is larger when the 56 Ni heating 
gets out, and therefore the temperature will be lower. 

2. The luminosity at early times is increasing faster than 
the radius expands, causing temperature to evolve bluer. 

3. The photospheric velocities are lower. 

4. The photospheric velocities evolve more slowly with 
time. 

2.7. The Minimum Explosion Time 

In additional to the qualitative correlations listed above, 
a single measurement of the photosphere velocity and color 
temperature during the rise can provide a strict, model inde- 
pendent, upper limit to the time of explosion before that mea- 
surement. This is seen by first rewriting equation (fT4l as 



L^47T,ja SB T c 4 /T c , 



(16) 



and then setting r c 
explosion began, to estimate that 



v c t e xp where f exp is the time since the 
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In general T c can be measured directly from the observations, 
but any velocity that can be measured will be from an ab- 
sorption feature shallower than the photosphere and generally 



Vph 



> 



Furthermore, t c will be greater than unity but also 



difficult to infer just from the observation. It is therefore use- 
ful to have a quantity that can simply be estimated directly in 
terms of observable quantities, 



1/2 



4ttv 2 a SB T c 4 



r ./2 ^J~W 



(18) 



Substituting typical values we find 



= 4.3 



10 42 erg s" 



10 4 K 



f / ph . V'days. (19) 
UO^ms- 1 / y V ; 



Therefore f m j n is an observable quantity that provides a model 
independent lower limit to the time of explosion^ Another 
reason that f m ; n is especially useful is that it only requires that 
the velocity and temperature of the SN be obtained at a single 
time. Thus when resources are limited, using equation ( fT9l 
is a helpful technique for learning a lot about the SN with 
minimal additional investment. 

This concludes our discussion of the qualitative features of 
rising SN I lightcurves. In general specific events may have 
details that we have not addressed, like non-spherical explo- 
sions or complicated velocity profiles. But by clearly spelling 
out the chain of reasoning behind our conclusions we pro- 
vide rules of thumb that can be used to build intuition, even in 
cases that are somewhat more complicated. To illustrate how 
our arguments can be used in practice, in the next section we 
apply them to a recently-discovered, well-studied SNe Ic. 

3. PTF 10VGV 

PTF lOvgv is a SN Ic that was discovered on 2010 Septem- 
ber 14.1446 (UTC time) with the Palomar Oschin Schmidt 4 8 
inch telescope (P48), by the PTF survey dCorsi et alj|2012l) . 
In a previous image taken on 2010 September 12.4830, it was 
not seem down to a limiting magnitude of R > 20.2. Follow- 
ing detection, the /?-band luminosity rises quickly to a peak 

10 days later. A single spectrum was taken « 2 days after 
detection. In light of what we have been discussing, it would 
have been ideal if the emission from shock heating could have 
been identified. It would have provided a measurement of the 
progenitor radius and the time of explosion, which could have 
been used to directly probe the 56 Ni deposition during the ris- 
ing lightcurve. Although it is unfortunate that this was not 
available, this event is ideal as a test case for exploring what 
can be learned when the time of explosion is not known. 

3.1. Radius Constraints 

Even though shock-heated cooling was not observed, PTF 
lOvgv did have an early detection of the rising lightcurve 
and upper limits in the time before this. Using this infor- 
mation, inter esting constraints on the progenitor's radius are 
still possible dCorsi et alj|2012l) . In Figure [3] we plot the ob- 
served /?-band absolute magnitude of PTF lOvgv (circles), 
along with upper limits (triangle). These all assume a distance 
modulus of 34.05 and galactic extinction of Ar = 0.445 mag 
dSchlegel et alj f 1998). We then calculate the shock-heated 
cooling according to equation ([TJ, which transitions to the 
plateau phase given by equation ((5]). We infer the corre- 
sponding /?-band absolute magnitude using equation (0, or a 
plateau temperature of ps 0.6 eV, along with P48-calibrated 
bolometric corrections provided by Eran O. Ofek (private 
communication). The resulting lightcurves are plotted along- 
side the data from PTF lOvgv in Figure [3] for a range of ex- 
plosion times and progenitor radii. The conclusion is that 

6 In AppendixlEl we derive the value of r c and show that it is not much 
larger than unity. For example, even for high temperature of 20, 000 K and 
Vph = 20,000 km/s at 10 days after the explosion 1 < r c < 6, implying f ra j n < 
fexp < 2.5r ra j n . For lower temperature and velocity the range is even smaller. 
Therefore t mm is not only a lower limit but also a rough approximation of f CX p- 
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FIG. 3. — The circles and triangles plot the data for PTF lOvgv from 
Corsi et al. 12012) for the detections and upper limits, respectively. The 
curves show theoretical calculations of the shock heating lightcurves using 
equations (TJ and 0, which are set to plateau at a temperature of ~ 0.6 eV, 
for a range of explosion times and radii as labeled. In all cases we fix 
E = 10 5 ' erg and M = 2Mq for the sake of comparison. The main conclu- 
sion is that without direct identification of the time of explosion, the radius 
upper limit constraint can vary by an order of magnitude or more. 

without knowing the time of explosion, the radius can only 
be constrained to be R* < 1-20^0. In general, the earlier the 
explosion time is, the more stringent the constraints on the 
radius. 

3.2. Bolometric Lightcurve 

Next we consider what can be learned about the mass and 
distribution of 56 Ni from the observed lightcurve. Before this 
can be done, there are two considerations that must be made. 
First, we need to convert the observed /?-band magnitudes to 
a bolometric luminosity. We solve for this using 



L^4t:^ct sb T c 4 /t c 



, yq(M r , & -M r -BC)/2.5 l 



(20) 



Since BC depends on T c , we can self-consistently find T c and 
thus L at any given time. For our estimates here we simply 
take t c ~ 1 and r c ~ r p h as given by equation ( fT3l to produce 
the thick, solid curves in Figure [4] Although not physically 
accurate, these estimates are adequate since the bolometric 
luminosity is found to be fairly robust. Once again, we must 
take into account that the time of explosion is not known, so 
we consider three example explosion times. 

The second consideration we must make before attempt- 
ing to associate the bolometric lightcurve with the 56 Ni dis- 
tribution is the impact of a diffusive tail. Assuming that the 
bolometric lightcurves in Figure [4] are representative of the 
56 Ni distribution, we can then ask, for a given luminosity at a 
given time what is the implied diffusive tail at earlier times? 
The dotted lines in Figure [4] show the diffusive tails originat- 
ing from various times using equation ( flOb . If the dotted line 
sits below the solid curve, this means that the diffusive tail is 
insufficient to explain the bolometric lightcurve at this depth, 
and thus the bolometric lightcurve represents direct heating 
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FIG. 4. — The bolometric luminosity of PTF lOvgv (thick, solid curves). 
From the top to bottom panel, we consider explosion times of 2, 3 and 5 days 
before the first detection. In all cases we fix E = 10 51 erg and M = 2M Q . For 
each lightcurve, at various times we consider what the diffusive tail coming 
from that depth would look like at earlier times (dotted lines) using equation 
fiOl From this comparison we conclude that if the SN occurred recently (top 
panel), then only the earliest heating may be explained by a diffusive tail and 
the majority of the rising lightcurve directly probes the 56 Ni distribution. On 
the other hand, if the explosion occurred further in the past (the middle and 
bottom panels), then the majority of the rising lightcurve can be explained 
as a diffusive tail due to deep deposits of 56 Ni. This example demonstrates 
the degeneracy between the time of explosion and the depth of 56 Ni if only a 
photometric lightcurve is available. 



from 56 Ni (as shown in the top diagram in Figure |2). On 
the other hand, if the dotted line sits above the solid curve, 
then the diffusive tail from that point overestimates the earlier 
bolometric lightcurve and luminosity at that time cannot be 
from direct heating of 56 Ni at the diffusion depth (as shown in 
the bottom diagram in Figure|2]). 

The conclusion from this comparison is that it is difficult 
to determine which of these three explosion times are more 
accurate from just this information. If the SN occurred re- 
cently (top panel of Figure IH, then only the earliest heating 
is explained by a diffusive tail and we would infer that the 
majority of the rising lightcurve directly probes the 56 Ni dis- 
tribution. On the other hand, if the explosion occurred further 
in the past (the middle and bottom panels of Figure [U, then 
the majority of the rising lightcurve is explained as a diffusive 
tail due to dee p dep osits of 56 Ni. This is consistent with our 
discussion in jj2.6l that there is a degeneracy between these 
two limits unless more information is available. 

Nevertheless, we can still try to constrain the mass and dis- 
tribution of 56 Ni as a function of the explosion time, the results 
of which are shown in Figure [5] In the top panel we plot the 
mass of 56 Ni inferred from equation (O. This is only plotted 
for depths at which the bolometric luminosity cannot be ex- 
plained by a diffusive tail using Figure [4] In the middle panel 
we plot the mass fraction X^, and in bottom panel the thermal 
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FIG. 5. — Inferred mass, mass fraction, and depth of 56 Ni in PTF lOvgv for 
different explosion times. The 56 Ni mass M56 is inferred from equation {6}, 
and we only plot M56 and X56 for times at which the bolometric luminos- 
ity cannot be explained as a diffusive tail according to Figure [4] The mass 
fraction is estimated as X56 pa M56 / AMjjff, where the diffusion depth is cal- 
culated according to equation {8). 
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FIG. 6. — Temperature and velocity evolution for PTF lOvgv when fitting 
a temperature of 6700 K and velocity of 16,000 km s at 2 days past first 
detection (shown as filled circles). In this particular case, the best fits occur 
for £51 =0.41M°' 72 . The curves are theoretical calculations assuming differ- 
ent explosion times as labeled, which demonstrates than an explosion time of 
~ 5 days before first detection is favored. 



diffusion depth using equation (0. 

3.3. Temperature and Velocity Constraints 

We next consider what additional information can be 
learned about PTF lOvgv from temperature and velocity mea- 
surements. U nfortunately, the av ailable data during the rise is 
rather limited. ICorsi et al.l (120 121) mention a 16, 000 km s" 1 Si 
II absorption feature (which is t ypically associated with the 
photosphere, Tanaka et al. 2008?) observed on 2010 Septem- 
ber 16, which corresponds to w 2 days after discovery. From 
the spectrum taken on that same day, by eye one can see that 
it roughly peaks at ~ 4300 A, which, assuming a blackbody 
spectrum, corresponds to a temperature of ~ 670 K. But we 
must b e cautious. In the study of SN 19941 by ISauer et al.l 
d2006l) . a spectrum with a similar peak is observed at 8 days 
post explosion, and their detailed models infer a temperature 
of ~ 10,000 K. As we describe next, if the correct tempera- 
ture were known, then tight constraints could be placed on the 
time of explosion and thus the other properties of PTF lOvgv. 
But given that the temperature cannot be extracted from the 
data without a detailed spectral modeling, which is beyond 
the scope of this paper, we separately consider both the cases 
of 6700 K and 10,000 K and discuss the implications. 

Using equation (fT9l we estimate the minimum time of the 
explosion. For T c = 6700 K, v p h = 16,000 km s -1 , and L « 
1.5 x 10 42 erg s" 1 (taken from Figure|4j, we find f m i n ~ 7 days. 
Therefore, if the cooler temperature is the correct one then 
the explosion must have occurred ~ 5 days or more before 
the first detection. On the other hand, for T c = 10,000 K we 
estimate f m j n ~ 3 days and the explosion must have occurred 
~ 1 day or more before the first detection. 

To test these conclusions, we perform a detailed fit of the 
temperature, velocity, and bolometric luminosity constraints 



simultaneously. This is done by varying E, M, and the time 
of explosion over a wide range of values and identifying 
which best fit the constraints. The result for either 6700 K 
or 10, 000 K is that a fit can only be obtained when E oc M° 72 
because this fixes the normalization of the velocity (see eq. 
lfl5l ). But for a given temperature only a very narrow range 
of explosion times are consistent with the data. In Figures[6] 
and [7] we plot the results of our fitting. These show that for 
a temperature of 6700 K or 10,000 K, the explosion much 
have occurred ~ 5 days or ~ 2 days prior to first detection, 
respectively. These match the values of t m \ n estimated before. 

Besides showing how the explosion time can be con- 
strained, Figures|6]and|7]highlight many of the general trends 
that we discussed in ^2]when considering how different explo- 
sion times impact various observables of the SN. For a more 
recent explosion time, the gradient in v P h is greater and the 
temperature is higher and decreasing. In contrast, for an ex- 
plosion time more distant in the past, T c is lower and increas- 
ing for a longer time. If merely a couple of data points were 
available, these trends could be identified in the observations, 
and even tighter constraints could be placed on PTF lOvgv. 
Nevertheless, it is powerful that even a single velocity and 
temperature measurement provides such stringent constraints. 
SNe la tend to have more early data available than SNe Ib/c, 
and we consider some of these events in forthcoming work. 

3.4. Progenitors 

Our conclusions from the previous sections show just how 
greatly our inferences about the progenitor can vary depend- 
ing on the explosion time. These are summarized as follows. 

1. If the temperature at 2 days past first detection is 
~ 6700 K, then the explosion occurred ~ 5 days or 
more before first detection. The radius is constrained 
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FIG. 7. — The same as Figure [6] but this time fitting a temperature of 
10,000 K. Now the best fits occur for E sl = Q3QMf 72 . In this hotter case, 
the explosion time is constrained to be ~ 2 days before first detection. 

to be < IRq and the 56 Ni must be located much deeper 
in the ejecta. Also the ejecta is likely of higher mass 
since it takes longer to reach the SN peak. 

2. If the temperature at 2 days past first detection 
is ~ 10,000 K, then the explosion occurred merely 
~ 2 days before first detection. The radius is con- 
strained to be < 6Rq and the 56 Ni must be located fairly 
shallowly. Also the ejecta is likely of lower mass. 

So we are left with two seemingly opposite conclusions about 
almost all the aspects of the ejecta just from opposite assump- 
tions about the time of explosion. We hope that a proper mod- 
eling of the available spectrum, which we do not carry out 
here, can provide the correct temperature and decide between 
these two options. We next discuss the implications for the 
progenitor star that produced PTF lOvgv in each case. 

Since the progeni tors of such SN Ic are hydrogen-stripped 
stars, the models of lYoon et al.l (120101) are a helpful starting 
place. The most striking feature of these progenitors is the 
inverse relationship between mass and radius that is shown 
in their Figure 12. If PTF lOvgv had a radius < IRq, it is 
strongly inconsistent with any of the binary models below 
Mf < 7M Q , where Mf is the final mass at the time of ex- 
plosion, which all have radii larger than ~ 2Rq. Their helium 
star models (i.e., Wolf-Rayet progenitors) have smaller radii 
and may be consistent with Mf > 6 - 8M Q , depending of the 
mass loss prescription employecQ. Other progenitors that are 
consi stent with such sm all radii are the gamma-ray burst mod- 
els of lYoon et all (120061) with Mf > 12M Q , but these probably 
are not applicable to PTF lOvgv. 

Combining the generally large mass inferred for a small- 
radius progenitor with an ejecta mass of ~ 3-5M Q (estimated 

7 In previous calculations of helium stars, such as by Woosley et al. 1 1995), 
smaller radii are generally seen. This differenc e is thought to be due to both 
different mass loss rates and updat ed opacities; [2?on et al. (201CJ) use opaci- 
ties from Islesias & Rogers ( 1996). 



from AMdiff at peak), the remnant mass must be > 2- 3M Q . 
This is near the range of t he maximum mass of normal NSs 
dLattimer & Prakashll2001l) . Would this imply that this event 
led to the formation of a black hole? Given the number of 
estimates that have been made throughout our analysis, it is 
difficult to make this conclusion with such certainty. But this 
does demonstrate how our analysis can make interesting con- 
nections between observations on one hand and the detailed 
simulations of progenitors on the other. 

For the case where R* < 6Rq, lYoon et aD (120101) show 
that such radii are naturally expected for progenitors in a 
mass range of Mf ~ 4-5M Q , and it is easier to reconcile 
the ejecta mass estimate with formation of a neutron star 
remnant. A less massive progenitor with a large radius has 
been predicted to have mixing of 56 Ni into the outer layers 
by Rayleigh-Taylor instability during the supernova explosion 
dHachisu et al.lll99TtlJoggerst et al.ll2009l) . similar to what we 
infer for a recent explosion (Fig. |5). So even in this case there 
are interesting correlations and implications for the explosion 
mechanism that can be investigated. 

Finally, it is worth discussing the inferences that can be 
made from the fact that PTF lOvgv was classified as a SN 
Ic. iDessart et al] (120121) studied the impact of 56 Ni mixing on 
the ejecta of SNe Ib/c, and found that it strongly impacts the 
observed spectral features. Chief among these are the He I 
lines, which require non-thermal electrons that are excited by 
7-ray lines fr om 56 Ni and 56 Co w hen the helium mass frac- 
tion is < 0.5 ( IDessart et al.112.01 lb . This means that the same 
ejecta with a significant amount of helium can produce either 
a SN lb if the 56 Ni maxing is strong or a SN Ic if the 56 Ni 
mixing is weak. If PTF lOvgv were constrained to have a re- 
cent explosion with strong 56 Ni mixing, we would be forced 
to conclude that its outer layers are very helium-poor. In the 
future it would be informative to repeat our analysis on a SN 
lb to see if shallow 56 Ni can be inferred to make the He I 
visible. 

4. CONCLUSIONS AND DISCUSSION 

We have discussed the photometric rising lightcurves of 
SNe I to investigate what can and cannot be learned from these 
observations. We provided a detailed summary of the vari- 
ous stages that determine the early lightcurve, and how their 
relative contributions are impacted by 56 Ni mixing and uncer- 
tainties in the explosion time. We then looked at a particular 
SN Ic event PTF lOvgv in some detail. Even with no direct 
constraint on the explosion time, and from just an /?-band pho- 
tometric lightcurve and a single estimate of the velocity and 
temperature of the ejecta, we argued that the explosion time 
could be constrained, which would have had a number of im- 
plications for the progenitor that we discussed. 

Our investigation demonstrates the kinds of connections 
that can be made between observations, explosion calcula- 
tions, and progenitor models. To facilitate such connections 
in the future, we summarize our conclusions below. We also 
highlight observational information that we deem especially 
important, and which should be considered of high priority, 
when investigating a rising lightcurve on limited resources. 

1 . A SN may exhibit a dark phase between the moment 
of explosion and the rise of the 56 Ni lightcurve. This 
means that extrapolating the 56 Ni lightcurve back in 
time is not a reliable method for estimating the time 
of explosion, and that without a known explosion time, 
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const raints on (such as for SN 201 lfe; lBloom et al.l 
I20T2I) are less stringent. 

2. Even though shock breakout may only emit at short 
wavelengths and may be too short-lived for detection in 
most circumstances, the UV/optical detection of shock- 
heated cooling phase can be just as useful for putting 
constraints on the progenitor radius. 

3. If caught when rising (eq. O), shock-heated cooling 
can also identify the time of explosion. Conversely, 
if the UV/optical rise is steeper than ~ ? L5 , then this 
argues that the shock-heated cooling is not being ob- 
served, and the explosion time is not well constrained. 
Due to a possible dark phase, data right before the 56 Ni 
lightcurve may not be sufficient to catch this rise that 
instead may occur ~ 1 — 6 days earlier. 

4. If the time of explosion is well constrained, then the 
56 Ni mixing can be extracted from the bolometric lumi- 
nosity in the following way. First use equation (O to 
find M^{t) assuming there is no contribution from the 
diffusive tail. Then at any time t ' find the diffusive tail 
contribution fromM 56 (f) at t < t' (eq. H0|). If this con- 
tribution falls below the actual luminosity thenMs6(f) is 
a reliable estimate (to within a factor of w 2) of the 56 Ni 
mass that mixed into AMdiff(f). Otherwise it is not. An 
upper limit on M^eit) can be derived then by finding the 
largest 56 Ni mass that produces a diffusive tail that falls 
below the observed luminosity. A useful rule of thumb 
is that if the luminosity rises faster than ~ f 2 , then it is 
dominated by the diffusive tail. 

5. If the time of explosion is unknown, having even a sin- 
gle temperature and velocity measurement during the 
rise can go a long way toward supplementing the pho- 
tometric data and provide strong constraints on the time 
of explosion using equation ( fT9] l. 



6. Having > 2 temperature and velocity measurements 
during the rise will allow the time evolution of each to 
be constrained in more detail. Knowing the time evo- 
lution of the temperature is especially useful because it 
provides a check of any conclusions about 56 Ni mixing. 

7. A future goal of observations should be to obtain sim- 
ilar coverage of SNe Ib/c with strong 56 Ni mixing. 
Such events would likely be cl assified as SNe lb u nless 
completely devoid of helium (iDessart et al.l 120121) and 
should never show a shock-heated cooling phase. This 
would be an important test of whether all the telltale 
signs of 56 Ni mixing are present as we have outlined. 

Although the scalings we have used are fairly robust, 
the numerical prefactors come from semi-analytic work 
that is only approximately correct. A useful future project 
would therefore be to calibrate our prefactors against more 
detailed numerical calculations, in particular the photo- 
spheric velocit y, color tempe r ature, and diffusion depth. As 
highlighted by IDessart et al.l d2012l) . the opacities for these 
ejecta can become complicated, so the various fits may only 
be applicable over different subsets of parameter space in 
composition and temperature. Nevertheless, the resulting 
collection of relations would be useful tools for quickly 
estimating properties of the explosions and the progenitors 
for future observations. 
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APPENDIX 



A. SHOCK-HEATED COOLING IN THE RAYLEIGH- JEANS TAIL 

At early times, when the lightcurve is dominated by shock-heated cooling of ejecta and before 56 Ni heating has become 
important, the lightcurve can rise in the optical/UV even if the bolometric luminosity is declining. This is possible as long as the 
optical/UV emission is in the Rayleigh- Jeans tail. In this limit 



^opt/UV °C 



1 



(Al) 



where r c is the color radius and r e is the optical depth at the color depth. Using the scalings from equations (20) and (21) of 
lNakar&Sarild20l0l) . we find that r c oc f° 83 and t c oc t 44 for a poly tropic index of n = 3. Combining this with T c from equation 
(O above, this gives L opt /uv oc f 1 5 as summarized in $2.21 Therefore when a lightcurve rises faster than f 1 5 it is strong evidence 
that something other than shock heating is provided the observed energy (although this conclusion is not foolproof, because in 
principle velocity gradients may differ from \> s oc p~° 19 as we assume). 

B. THE TOTAL CONTRIBUTION FROM THE DIFFUSIVE TAIL 

As mentioned in §!2.3l and $2. 41 the total luminosity from the diffusive tail in detail is not just from a single location. Instead it 
is integral over the diffusive tails from all depths. Such an integral can be approximated as 

l m- r\ (t^ dAMm c( ^ m ' lVit) dt' ran 
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where we take the limit of the integration to be the time of peak t pea ^ since this roughly corresponds to when the diffusion wave 
has made its way completely through the ejecta. Making use of the time-dependence of the diffusion depth from equation ([8]), 
the integral can be rewritten as 

In Figure [8] we plot the results of this integral for various distributions of X^. In all cases we use a constant luminosity from 
56 Ni heating (no exponential decay) for simplicity. In the first three plots we keep the shape of the increase of Xse the same, but 
vary its depth, increasing in depth from the top left to the top right to the bottom left. The top panel of each plot shows the various 
contributions to the lightcurve as labeled. The middle panels show the ratio of the tail contribution to the direct heating with a 
dotted line at a ratio of unity. In the bottom panel shows the distribution. In the bottom right plot we consider a shallower 
distribution of 56 Ni, as can be seen its the bottom panel. 

These example confirm our general conclusion that when the 56 Ni is more shallow or has a shallower slope, the impact of the 
diffusive tail is diminished. These examples also show that in the shallow outer layers exposed at early times, the diffusive tail is 
always at least comparable to the direct 56 Ni luminosity, so that cannot be inferred from the observations more accurately than 
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a factor of ~ 2. This uncertainty is not overly problematic because of the number of other estimates in this work and estimates of 
M56 within an order of magnitude will still facilitate our goals. Such a limitation should not be surprising, because the diffusive 
tail always leaks at least a little and contami nate s Lse- Our investigation of this effect also shows that the ratio Ltaii.totai /^56 
strongly depends on the prefactor in equation (IB2b . which comes from the power-law scaling of AM^ oc f. We generally find 
that this ratio is weaker with smaller r), demonstrating that L56 can be estimated more reliably in such circumstances. 

C. THE DIFFUSION DEPTH 

Here we d erive the depth of the diffus ion wave as a function of time, as was used in equation ((HJ). This roughly matches the 
derivation bv lRabinak & Waxmanl d201 ll) . with a few changes of numerical factors. 
The pre-explosion density profile of the star is approximated as 

P0(r ) = Pl/2^, (CI) 

where p\i 2 is the density at half the radius, 5 = 1 -ro/R* is the fractional radius, and n is the polytropic index. The fraction of 
ejecta mass lying above a radius ro is 



5 " = Ti f ' ^r 2 oPQ (r )dr Q « MjL S "^ , 
M ■/(«)«, «+l 



so that the mass depth is AM = S m M. 

The velocity profile from the passage of the SN shock is taken from M atzner & McKeel d!999t) . who use a n interpolation 
between the Sedov-von Neuman-Taylor and the GanderMan-Frank-Kamenetskii-Saku rii self-similar solutions dVon Neumann! 
1 1 9471: iSedovll 1 9591: iTavlol 1 950t iGandeF Man & Frank-Kamenetskiil 1 956t ISakuraill 1 9601) . Near the surface of the star this gives 



v *~~ a \m) {wJ s " (C3) 

where A v w 0.8, (3 = 0.19, and f p = pin/ p* is the ratio of the half-radius density to the average density = 3M/4irRj. In 
Appen dix ID1 we derive f p for some appropriate density profiles for the SN I progenitors we are interested in. iMatzner & McKeel 
d 19991) show that the final velocity of the fluid is given by v/(ro) = f v (ro)v s (ro), where f v m 2. 
The thermal diffusion length scale is 



D(S m ,t)tt x l " (C4) 
3np(d m ,t) 



where the density is given from continuity to be 



M fdv f \ l + l/n AM 
p - A^iXdS^) ~ /3 Airily ( } 

The diffusion depth is approximated as the location where D is roughly equal to the thickness of a given shell of material 

^W^T&T- (C6) 

do m \ + \/no m 



Setting Ar = D, we find 



Sn, 



2 / E- \ 1/2 / „ \jS/, f \ /3/a+l/")! < 1+1 /")/(l+l/n+/3) 



ft 3k M \MJ \3fJ U+l 



(C7) 



This expression is then used to derive equation ([8]). Note that AM^ = M5 m oc /^°' 04 for /3 = 0.19 and n = 3. In Appendix ID1 we 
find that f p ranges from rj 0.01 for SN Ib/c progenitors (this is what is used for the numerical factor in eq. O) to rj 1 for SN la 
progenitors, so the prefactor of equation ([8} varies by w 20% between these two very different scenarios. 

D. THE DENSITY STRUCTURE FACTOR 

Even though the diffusion depth is relatively insensitive to the d imensionless density struct ure factor f p , it is worth approximat- 
ing just how greatly it varies for different SN I progenitors. Using lMatzner & McKeel d 19991) . for a radiative envelope of constant 
opacity 

fpB JLJL(t^t] 4 ^a (Dl) 
Jp 144 1-/3 V k B J ' 

where 1-/3 = L^/L^a is the ratio of the stellar luminosity to the Eddington limit. For a Thomson opacity with hydrogen-deficit 
material this ratio is 

1-/3 = 0.51 (rir-)(^r-) ■ (D2) 



3Mr. 
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This then gives f p « 0.015 for a SN Ib/c progenitor. 
For degenerate electrons as in a WD progenitor for SNe la, the equation of state is P = Kp x+X l n . For the case of non-relativistic 

electrons n = 3/2 and K = 9.91 x I0 n p e 5 ^ 3 , and for relativistic electrons n = 3 and K= 1.23 x 10 15 /j,^ 3 , where \x e is the molecular 
weight per electron and K is in cgs units. In this case, 



/p = 



4nR 



GM 



(n+l)KR, 



3M 

where in the last expression we have assumed relativistic electrons and fj, e = 2. 

E. THE THERMALIZATION DEPTH 

Here we approximate the optical depth at the thermalization depth as a function of the observed color temperature, T c , the 
photospheric velocity and the time since explosion. The thermalization depth is defined as the outermost region which is in 
thermal equilibrium. In case that scattering dominates the opacity, and assuming that free-free dominates the photon generation 
(and opacity), this is the location where free-free emission generates within a diffusion time just enough photons to sustain thermal 
equilibrium (see N akar & Sarill20 1 01 for a detailed discussion). Namely, «bb(T c ) = fdiff"ph.ff(7c) where «bb(7c) ~ aT* /3k^T c is the 

blackbody photon density and n p h,ff ~ 3 x 10 36 s -1 cm~ 3 p 2 T c 1//2 x,- is the volumetric rate at which photons with temperature ~ T c 
are generated by free-free (in cgs units), a is the radiation constant, k% is Boltzmann's constant, p is the mass density at the 
thermalization depth and x-, is the average number of ionized electrons per atom. Now, the diffusion time is fjiff ~ (r c r c /c)(3/( \ + 
1 /«) while t c = Kpr c (3/(1 + 1 jri). Therefore, taking j3 = 0.19 and n = 3 and neglecting the difference between v c and v p h, which 
introduces completely negligible correction, we obtain: 



t c \ 7/ V v. \ v 7 * xl/3 



io 4 k7 V 104km / s / V da y 



(El) 



-1/3 2/3 

where we approximate x, =1 and k \ = 1 , which are reasonable values for the typical early temperatures of T c = 10,000- 
20,000 K. For lower T c both x, and k, (which are linear at each other) are smaller and so is r c . This estimate of r c assumes 
that bound-free absorption is negligible, otherwise r c is smaller. Finally, t c > 1 is always satisfied. If equation (IEU indicates 
otherwise, then its assumptions are not valid (e.g., scattering does not dominate the opacity) and r c « 1. 
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